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ABSTRACT 



• J3943 



A single linear vertical passive array is used in the 
•SOFAR' channel to determine the depth of a single under- 
water source at a constant range. The phase and amplitude 
weights applied to the array are determined by the linear 
minimum variance estimation technique. The resulting beam 
pattern is compared to the conventional rime domain beam- 
former. It was found that the linear minimum variance esti- 
mation technique of amplitude shading and phase weighting 
was significantly superior to the conventional beamformer. 
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I. INTRODUCTION 



There are many solutions to the problem of predicting 
the location of an underwater energy source. One common 
solution is the use cf a passive hydrophone which detects 
the pressure waves radiating from the source. The hydro- 
phone sensor is assumed to be omnidirectional and therefore 
incapable cf estimating direction. To provide direction- 
ality a series cf sensors are placed in a row to form a 
passive linear array. 

A familiar method of determining directionality is 
time-domain beamforming. In this principle, it is assumed 
that the source is far enough away so that the pressure wave 
appears to be a plane wave when viewed at the site of the 
receiving array. 

Thus a set of time delays are calculated for any direc- 
tion cf signal arrival, which, when applied to the receiver 
outputs causes them to be in phase and to reinforce when 
summed. -The resultant angular response to signals arriving 
from other than the nominated direction is then a function 
cf the array geometry, relative no the signal wavelength, 
and any weighting factors which have been applied to the 
receiver outputs. The effect is to generate a main 
receiving beam in the desired direction, with a series of 
undesired subsidiary sidelobes whose magnitude can be 
controlled to some extent by the choice of a suitable array 
geometry and the use cf amplitude weightings on the receiver 
outputs. 

In crder to determine location, three or more such 
arrays separated by a known amount may be used. 

This study is ccncerned with a single linear vertical 
passive array and the determination of the depth of a single 
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underwater source. The analysis is based on the following 
assumptions : 

• The underwater source is emitting continuously and at a 
monochromatic freguency. 

• Both the source and receiving hydrophones are stationary 
in space causing a constant range. 

• The range is sufficiently long so that the channel is 
filled with R-R (ref rac ted-refrac te d) rays. 

• There is no distortion introduced in the propagating 
medium so that the signals received at each sensor are 
identical except for constant delays. 

• The source signal and noise are independent and 
stationary gaussian random processes. 

• The speed of sound profile is triangular and symmetric 
with the deep sound channel axis at 1000 meters. The 
velocity gradient is -0.017 mete rs/meter /sec above 1000 
meters and +0.017 meters/meter/sec below 1000 meters. 
This profile gives a speed of sound at the surface of 
1500 meters/sec. 

• The speed cf sound profile is constant in the horizontal 
plane . 

• Only R-R rays are considered. All other rays have 
sufficient less that their effect is negligible. 

As opposed to conventional time-domain beamforming, this 
study makes no assumption cf planar wave fronts at the 
receiver site. Therefore the time delays applied to each 

receiver will not, in general, be a linear function of 
depth . 

Since it is desired to determine whether or not there is 
a source present at a specific depth the result will be a 
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will indicate sianal 



source 



binary decision. A "1" 

present; a "0" will indicate signal source not present. For 
the constant range there will be "N" test depths investi- 
gated for the signal source. The number of hydrophones in 
the vertical array will be "L". 

Fcr a single source at a given depth, the travel time is 
calculated from the depth to each hydrophone. This travel 
time is converted into a phase delay for each hydrophone so 
that after summation from all hydrophones a maximum output 
is achieved. This output is then passed through a squaring 
device, an integrator, and a threshold and flip flop device 
to give a "1" binary output. If the signal source is at a 
different depth and the same previous phase delays are used 
for each hydrophone then the output will be somewhat less 
than the previous maximum. The difference in depth required 
to achieve a '*0" binary output is the depth resolution of 
the system. 

The travel times for each hydrophone are calculated for 
each cf the "N" source depths to be considered. "N" will 
ordinarily be much greater than "L" so that the system will 
be overdetermined. An over det ermined system is one in which 
there are more equations than unknowns. The objective then 
is to calculate the phase angle and amplitude weight for 
each hydrophone so that a determination can be made indi- 
cating the presence or absence of a signal source at a given 
dept h. 

The method used to calculate the phase and amplitude 
weights is the linear minimum variance estimation technique, 
linear minimum variance estimators are optimum when compared 
with all ether estimators for gaussian problems. The method 
is directly applicable to overdetermined systems. 

The output of the summer is calculated using the linear 
minimum variance amplitude and phase angles assuming a 
source at one of the "N" source depths and no source at the 
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ethers. The calculation is repeated for each of the depths. 
The result, when plotted against source depth, will be 
referred to as the "beam pattern" of the array in this 
report. (Although similar, it should not be interpreted as 
the angular response of an array as in the conventional 
definition of a beam pattern. The conventional definition 
loses much of its utility when the wavefronts are not 
planar.) Ideally the beam pattern will be maximum at the 
desired depth and very small at ail other depths so that the 
binary "1" decision will be made for a source at the desired 
depth, and a "0" decision for sources at all others. This 
team pattern is compared with the depth beam pattern of the 
conventional time-domain beamformer mentioned above. The 
purpose of this thesis is to determine, as an initial inves- 
tigation, whether the linear minimum variance estimation 
technique, when applied to a linear vertical array, is 
useful in depth discrimination at long ranges in a ’SOFAR* 
type sound channel. 
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II. GENERAL THE ORY 

a. eai acoustics 

The propagation cf sound in an elastic medium can be 
described mathematically by solutions of the wave equation 
using the appropriate boundary and medium conditions for a 
particular problem. The wave equation relating the acoustic 
pressure 'p' to the coordinates , x', , y , , , z , r and the time 
't', may he written as 

,2 2 ,2 ,2 .2 

a p _ c ap + ap + ap 

TT , 2 2 2 (2.1) 

d t ax ay a z 

where ' c' is a quantity that has the general significance of 
sound velocity and may vary with the coordinates. 

One may approximate the solution of the wave equation 
using .ray theory: its body of results and conclusions is 

called ray acoustics. 

Officer [Ref. 1] describes the ray solution as a 
complete solution to any particular propagation problem 
within the validity cf the approximation of the Eikonal 
equation to the wave equation. For these approximations to 
be valid neither the amplitude of the wave nor the speed of 
sound can change appreciably in distances comparable to a 
wavelength. 

Thus the path of a ray through a medium in which the 
speed of sound varies with depth can be calculated by the 
application of Snell's law 



cos8/c = 1/c Q = a constant for any one ray 



( 2 . 2 ) 
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where ' e 1 is the angle or depression made with the 
horizontal at a depth where the speed of sound is 1 c*, and 
1 cq ' is the speed at a depth (real or extrapolated) where 
the ray would become horizontal. 

In a medium in which the velocity of sound changes 
linearly with depth the sound rays can be shewn to be arcs 
of circles, that is, to have a constant radius of curvature. 
Kinsler et al. [Ref. 2] give a simple and heuristic demons- 
tation of the circularity of rays in a medium with a linear 
sound speed gradient 'g'. The center of the circle which 
creates the arc lies at a depth where the sound speed extra- 
polates to zero. To understand this, consider a portion of 
a ray path with a local radius of curvature ’R', as illus- 
trated in Figure 2.1. Since the gradient ’g' for this case 
is 



g=Ac/AZ = (c 2 -c 1 ) / (d 2 -d = (c 2 -c 1 ) /R (cose 1 -cose 2 ) (2.3) 



where 1 a c* is the change in 
change in depth. It can be 
ture is given by 



sound speed and z' is the 
seen that the radius of curva- 



R = -c /g = ”C/(g ccs0) (2.4) 

The ray path is therefore a circle when • g • is constant 
because • R • is then constant. The center of curvature of a 
circle lies at the depth where 6 is 90 degrees, which corre- 
sponds to c=0. For the situation in Figure 2.1 the speed 
gradient is negative so that 'R' is positive. If the speed 
gradient were positive ' R' would be negative, and the path 
would curve upward. 

Once the radius cf curvature of each segment of a path 
is known the actual path can be traced graphically or 
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Figure 2.1 Circular Say Path 
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computed. If the initial angle of depression of a ray is 
01, then by referring to the geometry of Figure 2.1 the 
changes in both range and depth are 

A r= c x (sin (9 1) -sin (02) ) /(g cos(91)) (2.5) 



A z= c x (cos (02) -cos (01) ) / ( g cos (91)) (2.6) 

The sign convention for these equations is: downward, to the 
right, and depression angles below the horizontal axis 
positive. 

The symmetric triangular sound speed profile assumed in 
the introduction is similar to speed profiles encountered in 
the deep sound channel, sometimes called the SOFAR channel. 
The velocity minimum which occurs at the axis of the sound 
channel causes the sea to act like a kind of lens; above and 
below^the minimum, the velocity gradient continually bends 
the sound rays toward the depth of minimum velocity. A 
portion of the power radiated by a source in the deep sound 
channel accordingly remains within the channel and encoun- 
ters nc acoustic losses by reflection from the surface and 
bottom. These rays are called R-R (r efracted-ref racted) and 
since they have very low transmission loss, very long ranges 
can be obtained from a source of moderate acoustic power 
output. Thus an energy source at a specific depth will 
propagate energy in all directions but only the direction 
which is toward the receiving array and for which the rays 
are R-R is of interest. To determine the range of depres- 
sion angles which will yield R-R rays Snell’s law is used to 
give 
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9max = arc cos (c 1/ (cl- (g d) ) ) 



(2.7) 



where '©max' is in radians and 'cl* is the speed of sound a: 
the source depth 'd'. 

An example is given in Figure 2.2 of a single ray trace 
propagating in the SOFAR channel no show how equations. 2.5 
and 2.6 are used in determining range and depth. The ray 
path is broken up into arcs of circles as shown in the 
figure, and then by paying close attention to the previously 
defined sign conventions, the change in range and depth is 
found. Being more specific, for arc 1, 91 and 02 are both 

positive; for arc 2, 61 is positive and 92 is zero; for arc 
3, 01 is zero and e2 is negative; and for arc 4, 01 is 

negative and 92 is zero. By keeping a running total of all 

the depth and range changes it is possible to determine the 
total horizontal distance travelled and the depth at that 
distance. 

For the speed of sound profile assumed in the introduc- 
tion a computer generated ray plot is shown in Figure 2.3 
for a source depth cf 300 meters. As each ray propagates 
out from the source the triangular channel becomes filled 
with sound. .If a receiving hydrophone is placed a great 
distance away, a number of refracted propagation paths will 
exist, each having a different travel time and crossing the 
channel axis at different intervals. The path with the 
greatest excursion from the axis will have the shortest 
travel time. 

Officer [Ref. 1] shows that the travel time ' t' of a 
ray, which is an arc of a circle, is given by 



t 




de 

COS0 



( 2 . 8 ) 



and the travel time for each arc is 
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t 



1 log e 
g 



(2.9) 



tan( ( tt/4 ) + ( 9 2 /2 ) ) 



t an ( ( tT/ 4 T + ( 9 i / 2 ) ) 



Equation 2.9, applied using the same convention as equations 
2.5 and 2.6, determines the total travel time of a ray in 
the deep sound channel. 



B. AEBAY HCDEL 

The receiving linear vertical array is presumed to 
consist of 1 L* hydrophones as shown in Figure 2.4. It is 
assumed that the source is emitting energy at a constant 
frequency, ' f * , and amplitude, ’A', regardless of the depth. 
The source signal at the source is Aexp(j2 TT ft). The 
inherent received signal at the first hydropnone is 



xl(t) = A ' exp ( j2t f (t-t 1) ) (2.10) 

where *t1’ is the travel time from the energy source to the 
first hydrophone and ' A* 1 is the amplitude of the signal at 
the range cf the array. After passage through the amplitude 
weight * al 1 and a phase delay of ' ti», the signal on the 
first hydrophone at the input to the summer is 



yl(t)=A'alxl(t-Ti ) =A' alexp( j2irf(t-tl-xi)) (2.11) 

A time delay, for monochromatic signals corresponds to a 
phase shift 

9 = 2irft (2. 12) 

where *9' is the phase shift in radians and ’t 1 is the time 
delay in seconds. Thus equation 2.11 can be written as 
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y 1 (t) =A 'a lexp ( j(2 



(2. 13) 



where 1 = 2tt f h 1 * is zhe phase delay due zo zhe travel time 

from the source to the firsz hydrophone and '01* is the 
phase delay in the receiver on the first hydrophone. 

Combining all the hydrophones in the array in a summer 
gives as an expression for the array output 



Y(t) = l A'a k exp ( j(2irft-$ k -e k ) ) 
k — 1 



(2.14) 



where *<j> k ' represents the phase delay due zo the zravel time 
from the source to the ”k th” hydrophone and ' 6^ ' is the 
phase delay in the receiver on the ” k th” hydrophone. If 
the amplitude of the energy source is normalized by sezting 
A * = 1 , and equation 2.14 is written in terms of real and 
imaginary components, we have 



Y(t)= £a k cos(-<Ji k -e^)+jsin(-(j> v -9 v ) exp(j2Trft) 



k=l 



k k' 



(2. 15) 



When an energy source is at zhe ”q th” depth, we wish to 
have each receiving hydrcphcne's phase delay cancel out the 
effect of the zravel time from the source to it, such that 
* — 1> k~ € k * is squal to a multiple of This will put all 

signals into the summer in phase and thus maximize the 
signal gain for a source at zhe ”q th” depzh. From equation 
2.15 , 



£ cos ( -<J>kq“®kq^ = ^ ^ source present at q) (2.16) 

w ™ 1 



I ^ a ]<q ~^kq _9 kq^ = ^ (source present at q) 
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Note that -he firs- subscript on -he phase angle indicates 
the receiving hydrophene and the second subscript indicates 
the depth of the source. Thus ' j)j, ' would indicate the 

phase shift relating to the travel time from the "q th" test 
depth to the "k th" hydrophone in the receiving array. 

It is desirable for ’ Y (t) ' to be a minimum value for 
sources at other than the depth being investigated. Thus 
for each of the other * N-1 ' depths ' Y (t) ' is set to zero. 
This gives 'N-1' equations for the real terms of 'Y(t) ' set 
to ze ro 



J a kq cos (-fkm-ekq) =0 (source absent •> m/q) (2.18) 

k=l 

and 'N-1* equations for the imaginary terms of ’ Y ( t ) * set to 
zero 



L 

l a k sin tknm~ 0 :<q) =0 (source absent', m/q) 
k = l 4 4 



(2. 19) 



By using elementary tr igoncmetric identities, equation 
2.16 (real terms with source present at "q th" depth) 
becomes 



L 

J a kq [ cos (tkq) cos < 0 kq)“ sin (fkq) sin (9kq)J =1 



( 2 . 20 ) 



Equation 2.18 (real terms with source absent for each of the 
ether 'N-1' depths) becomes 






cosU km )COs( 9 kq ) - sinU km )sin(9 kq 



,)] =0 m=l , 2 , . . .N ; m/q 



( 2 . 21 ) 



Equation 2.17 (imaginary terms with source present at the "q 
th" depth) becomes 
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( 2 . 22 ) 



Y -a (sin(f )cos(6 )+ccs(J) ) sin (9 ) ) =0 

kq 1 kq kq 1 kq kq 

and equation 2.19 (imaginary terms with source absent for 
each of the other ' N-1* depths) becomes 



L r 

■kl 1 ^ ° ^ ^ 5 C ° S C 9 *q : 5 +C ° S ( ' *kn : ) 5 in ( 9 kq } ] = 0 



m=l,2,...N; mj!q 



(2.23) 
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the lowest (shallowest) source depth and equations increase 
in order after that until the last real and last imaginary 
equation correspond tc a source at the deepest depth. The 
resultant matrix equation becomes: 



o 

L 

o 




a iq cos9 : q 
a 2q COS0 2q 



a Lq COs9 Lq 



a lq Sin9 lq 

a 2q sin9 2q 



a T sin9 T 

. Lq L ^J (2. 24) 



Simplifying further, equation 2.24 becomes: 
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(2.25) 



where ' a,', *j_' , 'c_' , and represent the appropriate 

submatrices. 

Then, ty noting that the multipl icazicn of each element 
of a column by the same nonzero constant doesn't affect the 
solution, equation 2.25 becomes: 




(2. 26) 



Finally letting Z, A, and 9 represent the matrices in equa 
tion 2.26, we obtain 



_Z = A0 (2. 27) 

where a matrix is denoted ty a capitalized underlined 
letter. 

In summary, Z is the desired response of the vertical 
array to the 'N' source test depths. A represents known 
travel times from each of the 'N' source depths to each of 
the ' L' receiving hydrophones. jjj. is unknown. It is the 
phase and amplitude weighting which must be applied tc the 
vertical array in crder to realize Z_. ^contains '2L' 
unknowns. 
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Equaticr. 2.27 represents ’2N 1 equations. Since this 
system of equations is o ver det er mined (N>L) an exact solu- 
tion dees net exist. In order to make the best estimate of 
6 for the. desired response, the linear minimum variance 
estimation technique is used. 

C. LINEAB BIN I HUM 7 ABI&NCE HETHOD 

Equation 2.27 represents a noise free environment. If 
noise were present it would become 

_Z = A0+n_ (2. 28) 

with 'n/ a ^ 2N ,, element column vector. This represents the 
noise at each source depth. *_Z_* is a linear function of 
and is called the oberserva ticn. is a ”2Nx2L" modula- 

tion or observation matrix which is known, and ’e.' is the 
"2L" element random parameter vector which is to be esti- 
mated. Assume the first and second moments of • < 9* and 'n_' 
are qiven by 



E(0) = He V ar (8) = V 0 

and 



(2.29) 



E(n) =0 Var(n) = V n (2.30) 

where ' E* represents expectation or first moment and 'Var* 
represents the covariance matrix. It is assumed that the 
parameter *6' and the noise 'n.' are uncorrelated. 

A restriction imposed is that the estimate must be a 
weighted linear combination of the observations: 
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(2.31) 



§ = b + £.Z 

Li 

where ,A| indicates estimate. The objective is tc select 
•M and 'B_' in order tc minimize the error variance. Such 
an estimator is the linear minimum variance estimator; it is 
the test, in the sense of minimum-error-variance linear 
estimators. 

Another restriction is that the estimator be unbiased; 
in otherwords it is required that the expected value of the 

A 

estimator '0 ' is equal to the expected value of the param- 
Li 

etar *0'. Thus, 



E(9 l ) = fc + BE(z) =E(0) = jj 



(2.32) 



yielding 



b=b 0 -EAy 0 



(2.33) 



Substituting this result in equation 2.31 gives for the 
unbiased linear estiiiator 






(2. 30) 



Note that since the estimator is unbiased, the estimation 

A 

error '0 =6-6 ' is zero mean. The next step is tc select 

-E L 

'_B' in order to minimize the error variance. However, this 
optimization problem is ill-defined because the error vari- 
ance is a matrix. Therefore in order to introduce a scalar 
goodness measure the sum of the variances of each component 
of '0,' is minimized. This is the sum of the main diagonal 
terms of the covariance matrix and is defined as the trace 
of the matrix. 
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(2.35) 



2N 



tr (Var(0^ 
— L 

where ' tr ' 
minimize the 



)) = l V ar((9„) 
= 1 _t " 



indicates trace. 'B_* is 
trace of the error variance. 



then 

cr 



selected to 



min tr (Var(0J) = min tr (E (9 ei) ) (2.36) 

B L B L L 

T 

where ' ' indicates the transpose. The following problem is 

then obtained by substituting equation 2.34 into equation 
2.36: 



min tr ( Var(8_) ) =min tr (E( (9 -u g -3( Z-Au g ) ( 9-u 3 3( Z-Au Q ) T ) ) (2.37) 

B 3 

It is well known [Ref. 3] that equation 2.37 is minimized 
when 



Cov (9 ,Z)-BVar (Z) =0 






(2. 38) 


where ' Ccv(e,Z) ' is 


the covariance 


matrix of the 


u n kn ow n 


parameters and the 


cbserv aticns . 


Denoting the 


optimum 


filter by *_B ** , then 


if 






B*=Cov (6 , Z) (7ar (Z) ) 


-1 




(2. 39) 



a minimum is achieved for the sum of the squares of the 
errors. 

Osing equation 2.28 for , < Z_* , the covariance of '9.' and 
•Z* becomes 



Co V (9, Z) =Cov (G,A9 + n)=V g A T 



(2.40) 



since * 0 * and • n* are un cor related. The variance of • Z' is 
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(2. 41) 



Var (Z) = Var (A0+n) =A.I 0 -i* + I n 
Substituting these into equation 2.39 gives 

B* = V Q A T (AV Q A T +V n )' 1 (2.42) 

and the linear minimum variance estimator is 



£lMV = (2.43) 

By utilizing a matrix inversion lemma [Bef. 3] equation 
2.43 becomes 



— LMV 






(2. 44) 



The advantage of this equation over equation 2.43 is the 
size of the matrix to be inverted. In equation 2.43 the 

matrix has dimensionality '2N' while in equation 2.44 its 
dimenionality is only ’ 2L* . Thus the advantages of the 
linear variance estimator are the ease with which they are 
derived, the mathematical tiactabilizy of the linear form, 
and the minimum amount of stochastic information required 
for development. An interesting characteris tic is that the 
linear minimum variance estimate is the orthogonal projec- 
tion cf ’©* onto the space spanned by the observation ' Z* . 
Eecause of these factors this estimator is a popular form 
for estimating unknowns in cverde termined equations. 

For this thesis it is assumed that the noise samples are 
uncorrelatea and identically distributed so that: 



V = 



n 



2 

a I 



(2.45) 
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No previous knowledge is assumed about 'S'. This implies an 
infinite variance matrix which is represented as: 



Ie 1= o 



a nd 



if 0 =i 



(2. 46) 



The linear minimum variance estimate given by equation 2.44 
is then 



i L „f <i T i> ‘Vi 



(2.47) 



By determining ' G ' the phase and amplitude weights are 

“LMV 

found for a signal source on the 'q zh' depth. Recall that 



— LMV 



— ‘ “ 




Mi M ’ 


h 




a lq COS9 lq 


9 2 




a 2q COs9 2q 


9 l 


= 


a La cos9 Lq 


9 l+i 




a lq 31n6 lq 


9 2L 




a. sinQ- 

-q w 




L. -i 



(2. 48) 



Upon solving, this equation gives for the phase delay 



9 =arctan(9„ /9 ) where m=1 to L 
mo L+m m 



(2.49) 



and amplitude weigh- 



a =9 /cos (9 ) where m=1 zo L 

mq m mq 



(2.50) 



When these amplitude weights and phase delays are applied to 
the vertical linear array a resulting beam pattern is formed 
which in the absence of noise is: 
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( 2 . 51 ) 



Z = A 0 



The resulting beam pattern 'Z_ ' can then be compared with 

the desired beam pattern ' Z_' , as well as with a conventional 
beam pattern ' Z 1 * obtained using linear phase shifts across 
the array aperture selected to "steer" the array to the 
dominant arrival angle for the selected source depth. 
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Figure 2.2 Single Ray Path Plot In Triangular SOPAR Channel 
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Figure 2.3 Ray Plot For Assumed Sound Channel 



Source Depths Vertical Linear Array 
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Figure 2.4 array Model 



III. EXPER I MENT AL PROCEDURE 



A. EASIC ASSUMPTIONS 

The speed of sound profile given in the Introduction was 
used to test the techiniaue, with 'L', the number of source 
depths chosen as 20. The source depths are at 220 meters 
and every 20 meters thereafter to 600 meters inclusive. 
'N ' , the number of hydrophones used rn the vertical array 
was chosen as 5. The hydrophones are placed at depths of 
100, 200, 200, 400, and 500 meters. Thus, there are 40 

equations, two for each source depth, and 10 unknowns, two 
for each hydrophene. A range of 200 kilometers (km) assures 
that the deep sound channel is filled with sound over the 
aperture of the vertical linear array. A frequency of 100 
Hz provides good resolution in the beam pattern without 
introducing alias mainlobes at the selected range across the 
depths of investigation. 

E. »A/ MATRIX CALCULATION 

Since the gradient of the sound velocity, ’g' is 
constant in the area of the source depths, and the speed of 
sound at the surface is known, equation 2.3 is used to solve 
for '02*, the speed of sound at the source depth. Equation 
2.7 is used for each source depth to calculate the maximum 
initial depression angle which yields R-R rays. In deter- 
mining travel time, only depression angles from each source 
which are positive (downward) and yield R-R rays are used. 

Beginning with the first source depth of 220 meters, an 
initial depression angle of 0.0 degrees is selected. The 
ray path is then calculated using equations 2.5 and 2.6 and 
broken into a series of arcs as in figure 2.2. Equation 2.9 
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is used to determine the travel time for each arc in the 
same manner. Each arc's horizontal range, travel time, and 
depth are summed. When the summed horizontal range reaches 
200 km the summation process ends. The travel tim- and 
depth of the ray at this horizontal range is then known. 
The same procedure is repeated for an initial depression 
angle of 0.1 degrees and every increment of 0.1 degrees 
thereafter until the maximum depression from equation 2.7 is 
reached. The final ray path is at this maximum depression 
angle . 

The same procedure is repeated for each of the other 19 
source depths. 

Thus, for each initial angle from each source depth 
there is a ray which has a travel time and a depth when it 
reaches the horizontal range of 200 km. Since it is the 
profile of the sound pressure wave which impinges on the 
vertical array which is of importance, a constant can be 
subtracted from these calculated travel times. This 
constant is selected to be the travel time for the source 
depth of 220 meters which has an initial depression angle of 
0 degrees. It is subtracted from each of the travel times 
making the resultant travel times relative with respect to 
the ray which has a 0 degree depression angle from the 220 
meter depth. The program and its listing which calculates 
the relative travel times and the depths of these rays at 

the horizontal range of 200 km is given in Appendix A. 

A plot cf the relative travel times versus depth for the 
220 meter source depth is shown in Figure 3.1. Figure 3.2 

displays the plot for the 380 meter source depth. Negative 

relative travel times in the plots indicate that the overall 
travel time is less than the reference. These rays arrive 
at the 200 km horizontal distance before the reference ray. 

Since the receiving hydrophones are at set vertical 
positions (100 , 200, 300, 400, and 500 meters), an 
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interpolation is done to determine relative travel times to 
them from each source depth . The interpolation program and 
its listing is in Appendix B. Sometimes more than one R-R 
ray travels from the source depth to a hydrophone. when 
this occurs, the ray which arrives first is used in the 
calculation of relative travel time to that hydrophone. 

Eguation 2.12 determines the phase shift relating to the 
relative travel times. The 'JJ matrix is formed by taking 
the appropriate sine and cosine values as in eguation 2.24. 
The ' J_* matrix is *4 0 by 10*. 

C. »Z/ 1ATBIX 

Referring to equation 2.26, the matrix is a *40 by 

1' column vector. It is the desired beam pattern. The 
bottom 20 rows give the imaginary terms and are set to zero. 
The top 20 rows represent the value of the real terms at 
each source depth. Therefore each of the top 20 rows is set 
to zero except for the row containing the source. It is set 
to 1 . For example, if the source is at 220 meters then only 
the top row is set tc 1 . If the source is at 380 meters 
then only the ninth row is set to 1. 

D. BESOITIHG BE AH PATTEBH 

1 . Osinq the Lin ear Mi nimum Variance Method 

/v 

• e T ' is calculated using equation 2.47. The 
resulting beam pattern , Z > * is calculated using equation 
2.51. The program which calculates the ' A 1 matrix, uses it 

/N ^ 

in determining 1 aa< ^ then calculates *_Z_' is given in 

Appendix C. The program listing is also included. 
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2 • Osin 2 L ine ar Phase S hift s 



The conventional beam pattern is determined by using 

A 

equation 2.51 where is calculated by approximating the 

plot of relative travel time vs. depth by a straight line at 
the receiving hydrophone depths. For example Figure 3.3 
represents this plot for the 380 meter source depth. The 
straight line is determined by a least squares linear 
regression which mirimizes the sum of the squares of the 
deviations of the actual data points from the straight line 
cf best fit. Note that only data points which are on the 
dominant curve are used in calculating the straight line. 
From the straight line, relative travel times to the 
receiving hydrophones are calculated. The relative travel 
times for the 380 meter source depth are given in Table I. 
They correspond to a plane wave arrival angle of 3.73 

A 

degrees. is determined by converting these relative 

travel times to phase delays using equation 2.12 and then 
taking the appropriate sine and cosine values of these phase 
delays as in equation 2.24. The amplitude weights are 
initially assumed to be unity. 

A second method for obtaining the conventional beam 
pattern is calculated by the same procedure except the 
amplitude weights which are determined by equation 2.50, the 

/N 

• amplitude amplitude weights. 



’© 

-LMV 
hydrophone. 



are applied to each 
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TABLB 


I 


Relative Travel Tiaes For 


380 Heter Source Depth 


Hydrophone Depth 


Relative Tr aval Time 


100 meters 


-0.07565509 sec. 


200 meters 


-0.07131599 sec. 


200 meters 


-0.06697690 sec. 


400 meters 


-0.06263780 sec. 


500 meters 

i 


-0.05829871 sec. 
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Figure 3.1 Relative Travel Time vs. Depth (220 meter sourc 
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Figure 3.2 Relative Travel Time vs. Depth (380 meter source) 
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Figure 3.3 Straight Line Approx, of Rel. Trav 



IV. RESULTS 



&. EXACT SOLUTION 

An exact solution is derived for the bean pattern if the 
number of receiving hydrophones equal the number of source 
depths. For example, when the 5 receivers are used to 
discriminate between 5 source depths (220, 240, 260, 280, 

and 300 meters) there are 10 equations with 10 unknowns. 
Figure 4.1 is a plot of the resulting beam pattern with the 
source at the shallowest depth. Note that because of roun- 
doff errors in the 'IMSL' subroutines there is a small value 
for the resulting beam pattern at the non-energy source 
depths under investigation. 



B. FOUR DEPTHS NITH TWO RECEIVERS 



Fcr source depths at 220, 
receiving hydrophones at 100 
equations with 4 unknowns, 
resulting beam pattern with a 
Figure 4.3 is the plot for 
depth . 



280, 300, and 320 meters with 
and 200 meters, there are 8 
Figure 4.2 is a plot of the 
source at the 220 meter depth, 
the source at the 280 meter 



C. TWENTY DEPTHS WITH TWO RECEIVERS 

For all 20 source depths with receiving hydrophones at 
100 and 200 meters, there are 40 equations with 4 unknowns. 
Figure 4.4 is a plot of the resulting beam pattern for a 
source at the 220 meter depth. Figure 4.5 is a plot for the 
source at the 380 meter depth. 

A 

In order to determine if the •9' calculated in this case 
is the best an alternate method is devised. Four source 
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depths (220, 240, 260 meters, and another source test depth) 

with the original source at the 220 meter depth and 

receivers at 100 and 200 meters are used. The beam pattern 
v 

is calculated each time with a different source test depth 
substituted for the fourth source depth. The 5 best 
resulting beam patterns are selected along with the 8 source 
depths (220, 240, 260 meters, and the 5 test depths which 

created the 5 best beam patterns) . Then, using these 8 

A 

source depths, '9/ is determined for the two receivers. 

A 

This * 9 1 is applied to the two receivers and the beam 
pattern obtained for all 20 source depths. 

The resulting beam pattern obtained by this alternate 
method isn't as good as the beam pattern obtained by using 

A. 

all 20 source depths in the determination of '9'. 

D. TWENTY DEPTHS WITH FITE RECEIVERS 

For all 20 source depths with all 5 receiving hydro- 
phones there are 40 equations with 10 unknowns. Figure 4.6 
is a plot of the resulting beam pattern with the source at 
the 220 meter depth. Figures 4.7, 4.8, and 4.9 are the 

plots for the source at the 360, 380, and 400 meter depths 

respectively. 

E. CCNVEBTIOHAL BEA 8FORMEB 

Figure 4.10 is a plot of the beam pattern for a conven- 
tional beamformer using linear phase shifts across the array 
with the source at the 380 meter depth and the amplitude 
weights set to unity. All 20 source depths and 5 receiving 

hydrophones are used. Note that in figure 4.10 that -‘-here 

is less than 1 db discrimation between each of the source 
depths. Figure 4.11 is the plot obtained for the amplitude 
weights set to values determined by aquation 2.50. 



41 



F 



RANGE OF 250 KILOMETERS 



The calculations were repeated for a range of 250 km 
using the same 20 source depths and 5 receiving hydrophones. 
Figures 4.12 and 4.13 represent the plots of relative travel 
times versus depth fcr the 220 and 380 meter source depths 
respectively. Figure 4.14 represents the straight line 
approximation of the relative travel times for the 380 meter 
depth. Note that in this figure the relative travel times 
are represented by two straight lines; the upper line repre- 
sents linear phase shifts of the slower travel times for the 
conventional beamformer while the lower line represents 
linear phase shifts of the faster travel times. Tables II 
and III are the straight line interpolations of these slower 
and faster travel times which correspond to arrival angles 
of 4.04 and -3.16 degrees respectively. Figures 4.15 and 
4.16 are the resulting beam patterns for the conventional 
beamfcrmer for the slower travel times using unity amplitude 
weights and linear minimum variance amplitude weights 
respectively. Figures 4.17 and 4.18 are the beam patterns 
for the faster travel times. 

Figures 4.19, 4.20, 4.21, and 4.22 represent plots of 

the beam pattern for the source at the 220, 340, 360, and 

380 meter depths respectively. 
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T ABLE II 






Slower Bay Travel Times For 

Source 


250 km Range 


and 


380 m 


Hy dro phone 


Cepth 


Tr, 


avel 


Time 


100 meters 




-0 . 09949109 


sec . 




200 meters 




-0. 09478615 


sec • 




300 meters 




-0.09008 121 


sec . 




400 meters 




-0.08537627 


sec • 




500 meters 

1 




-0.08067133 


sec . 


J 












1 




TABLE III 






Faster Travel Times For 250 km 


Range and 380 m 


Source 


Hydrophene 


Depth 


Relative Tr 


avel 


Time 


100 meters 




-0. 10 100745 


sec. 




200 meters 




-0. 10468763 


sec . 




300 meters 




-0. 10836781 


sec • 




400 meters 




-0. 11204799 


sec . 




500 meters 

i 




-0 . 11572817 


sec . 


i 
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Figure 4.1 Beam Amplitude vs. Source Depth (5 depths 5 receivers) 
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Figure 4.2 Beaa Pattern (4 depths 2 receivers, source at 220 m) 
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Figure 4.4 Beam Pattern (20 depths 2 receivers, source at 220 m) 
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Figure 4.5 Beau Pattern (20 depths 2 receivers, source at 380 m) 
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Figure 4.6 Beam Pattern (20 depths 5 receivers, sourc 
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figure 4.7 Beam Pattern (20 depths 5 receivers, source at 360 m) 
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Figure 4.8 Beaa Pattern (20 depths 5 receivers, source at 380 m) 
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Figure 4.9 Bean Pattern (20 depths 5 receivers, source at 400 m) 
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Figure 4.11 Conventional Beam Pattern (non-unity amp. vts. sourc 
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Figure 4.12 Relative Travel Time vs. Depth (range=250 km, source at 220 n) 
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Figure 4.13 Relative Travel Time vs. Depth (range=250 km, source at 380 m) 
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Figure 4.15 Conventional Beam Pattern (R=250 km, d-380 m, slower times 



oro 




ao 


to 




CN 


o 


CO 


CO 




CN 


o 












o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


d 


d 


o 


o 


o 



< tsi 



59 



igure 4*16 Conventional Beam Pattern (H=250 km, d-380 m, slower times, LMV a) 
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Figure 4.18 Conventional Beam Pattern (B=250 km, d-380 m, faster times, LMV a) 
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Figure 4.19 Beam Pattern (R=250 km, 20 depths, 5 receivers, source at 220 m) 
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0.24-1 
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Figure 4.21 Beam Pattern (R=250 km, 20 depths, 5 receivers, source at 360 m) 
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Figure 4.22 Beam Pattern (R=250 km, 20 depths, 5 receivers, source at 380 m) 



V. DISCUSSION OF BESDLTS AND RECOMMENDATIONS 



For many of the cases examined, the linear minimum vari- 
ance estimation technique gives high resolution in the 
"depth beam pattern" for sources at long range. Figures 4.2 
and 4.3 represent beam patterns for four depths with two 
receivers where the main lobe is at the desired source depth 
with a beamwidth of 30 meters. For 20 depths with 2 
receivers the system is highly overdetermined and the beam 
pattern portrayed in figure 4.4 has its main lobe 80 meters 
from the source depth. At the source depth the strength of 
the beam pattern is 3.6 db down. For figure 4.5 the main 
lobe is on the source depth with a beamwidth of 250 meters. 

For all 20 source depths with 5 hydrophone receivers the 
results range from a beam pattern having its main lobe on 
the source depth (figure 4.7) with a beamwidth of 20 meters 
and a secondary lobe 12.1 db down to a beam pattern having 
its main lobe 40 meters from the source depth (figure 4.22). 
For the team pattern in figure 4.22 the width of the main 
lobe is 80 meters and the strength of the beam pattern at 
the source depth is 0.8 db down. 

When using the conventional beamformer, the beam pattern 
results with amplitude weights determined by the linear 
minimum variance estimation method were superior to the beam 
pattern determined by using unity amplitude weights. 
However in all cases the conventional beamformer was signif- 
icantly inferior to the beam pattern determined by the 
linear minimum variance estimation technique in terms of 
depth resolution. 

The linear minimum variance estimation technique could 
not be used for all ranges with the chosen sound speed 
profile. For ranges of 187 km and 235 km the relative 
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travel times produce an *kj matrix which is so ill- 

conditioned that the •IMSL* subroutine ' LINV2F' cannot 

. T 

determine an accurate inverse of 'A^ A'. 

There is sufficient evidence from this initial investi- 
gation that the linear minimum variance estimation technique 
applied to a long linear vertical array can yield high reso- 
lution depth information about passive sources at vary long 
ranges. However, further investigation is needed before any 
field tests are in order. Recommendations for further 
research are: 

1. The use of a more realistic speed of sound profile, 
preferably one actually characteristic of the ’ S CFAR* 
channel . 

2. The use of mere hydrophones in the vertical array. 
As more receivers are used the beam pattern should 
approach the desired beam pattern more closely (a 
less over determined system should yield smaller total 
minimum mean square error). 

3. Investigation into causes for the peak response 
falling at other than the desired depth and techni- 
ques for correcting this problem when using the 
linear minimum variance estimation techinique. 

4. Alternate assumption for choosing the ray paths to 
include in the *_A_' matrix. For example one might 
choose the rays with the greatest intensity instead 
of the shortest travel time. 

5. The possibility of estimating range as well as the 
depth of a passive source with linear minimum vari- 
ance estimation techniques. 

In conclusion, the linear minimum variance estimation 
technique of beamforming was signif icantl y superior to the 
conventional beamformer. High resolution depth information 
about passive sources at long ranges is provided. 
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APPENDIX A 



RELATIVE TRAVEL TIHE CALCULATION 

This program calculates the relative travel times from 
each source depth to the horizontal range and the depth of 
the ray at this range. To save paper, only the two shal- 
lowest source depths are listed. 
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APPENDIX 3 

INTEBPOIATICH OP RELATI7E TRA7EL TINE CALCOLATIONS 



This program does an interpolation 
generated in Appendix A. The relative 
interpolated for the receiving hydrophone 



of the 
tr avei 
depths . 



output data 
times are 
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RECEIVING HYDROPHONE DEPTHS ARE SELECTED. IMSL ROUTINE 'ICSEVU* 
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APPENDIX C 

RESULTING BEAM PATTERN POR CALCULATED WEIGHTS 

This program uses the output of Appendix B to calculate 
the 'A_' matrix. From this, the amplitude and phase weights 
are determined by the linear minumum variance estimation 
technique. These weights are applied to the array and the 
beam pattern is obtained. 
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